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Abstract 



Consider the setting of randomly weighted graphs, namely, graphs whose edge weights are 
chosen independently according to probability distributions with finite support over the non- 
negative reals. Under this setting, properties of weighted graphs typically become random 
variables and we are interested in computing their statistical features. Unfortunately, this turns 
out to be computationally hard for some properties albeit the problem of computing them in the 
traditional setting of algorithmic graph theory is tractable. For example, there are well known 
efficient algorithms that compute the diameter of a given weighted graph, yet, computing the 
expected diameter of a given randomly weighted graph is #P-hard even if the edge weights are 
identically distributed. 

In this paper, we define a family of properties of weighted graphs and show that for each 
property in this family, the problem of computing the fc*'' moment (and in particular, the 
expected value) of the corresponding random variable in a given randomly weighted graph G 
admits a fully polynomial time randomized approximation scheme (FPRAS) for every fixed k. 
This family includes fundamental properties of weighted graphs such as the diameter of G, the 
radius of G (with respect to any designated vertex) and the weight of a minimum spanning tree 
of G. 
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1 Introduction 



In the traditional setting of graph algorithms, the input is typically a graph G whose edges are 
often associated with some weights. In most applications, G represents a real-life network and the 
edge weights correspond to some attributes of the network's links which are assumed to be known 
when one wishes to apply some graph algorithm to G. Unfortunately, in many scenarios these 
attributes of the real-life network cannot be determined. Still, however, it is often believed that 
the attributes of the network's links are governed by some known probability distributions. For 
example, the latency along each communication link in the Internet backbone is usually assumed 
to be a random variable, rather than a fixed value that can be determined a priori. The same 
situation appears in the emerging field of networking called delay tolerant networks (DTNs) [13] 
that includes sparse ah-hoc networks |35^ S] , space exploration networks [9j , submarine networks 
|31j . and sensor networks. 

Our goal in this paper is to advance the theoretic foundations of graph algorithms operating 
in the context of edge weights that obey some specified probability distributions. In this context, 
various properties of the graph become random variables and we wish to design better algorithms 
for the problems of computing the statistical features of these random variables. This turns out to 
be a non-trivial task even for basic and fundamental graph properties which are easy to compute 
in the traditional (deterministic) setting. Similar challenges were previously tackled in many works 
(see the survey of Ball et al. [3J), however the computational angle of these problems received little 
treatment. 

The model. A randomly weighted (RW) graph is a graph0 G in which the edge weights are 
independent random variables with finite support over the non-negative reals. Specifically, every 
edge e G E{G) is associated with some positive integer m(e) and with some non-negative reals 
Wl, . . . , We"^^^^ and p\, ■ ■ ■ ,p^^^\ where X^^fi^ p\ = 1, such that the weight of e is w{e) = Wl with 
probability p\ independently of all other edges. The reals W}, . . . , W^^'^^ are called the phases of 
edge e. 

A special subclass of RW graphs that plays a major role in this paper is that of biased coin 
weighted (BCW) graphs in which m(e) = 2 for all edges e G E{G). It will be convenient to slightly 
change the notation for BCW graphs: edge e G E{G) is said to take on its heavy phase with 
probability p^ and on its light phase with probability 1 — Pe, where > > 0. A BCW 
graph G in which p^ = p, = , and = for all edges e G E{G) is referred to as an 
identically distributed weighted graph. 

^ Unless stated otherwise, all graphs in this paper are assumed to be finite and undirected, though not necessarily 
simple. The vertex set and edge set of graph G are denoted by V{G) and E(G), respectively. 
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Weighted graph properties. Let Q be the collection of all (not necessarily simple) connecteqj 
graphs G with non- negative edge weights w : E{G) — t- M>o. Throughout, we think of a weighted 
graph property as a function X that assigns a non-negative real X{G) to each graph G ^ Q. A 
weighted graph property X : Q ^ M>o is said to be distance- cumulative if it satisfies the following 
requirements for every graph G ^ Q: 

(Rl) If G' E ^ is the graph obtained from G by multiplying all edge weights by some factor 
r G M>o, then X{G') = r ■ X{G). 

(R2) li G' ^ Q is the graph obtained from G by increasing the weight of some edge e G E{G) by 
an additive term r G R>o, then X{G) < X{G') < X{G) + r. 

(R3) If w{e) =0 for some edge e G E{G), then X{G/e) = X{G), where G/e £ Q is the graph 
obtained from G by contracting the edge e. 

(R4) If e, e' G are parallel edges and w{e) > w{e'), then X{G - e) = X{G), where G-eeG 

is the graph obtained from G by removing the edge e. 

(R5) If Af(G) is strictly positive, then X{G) > mm{w{e) | e G -E;(G)}. 

(R6) X{G) = if and only if the graph obtained from G by removing all (strictly) positive weight 
edges is connected. 

Requirements (R1)-(R5) and the sufficiency direction of requirement (R6) are naturally sat- 
isfied by weighted graph properties in which the edge weights correspond to time delays, routing 
costs, etc. The necessity direction of requirement (R6) is slightly more restrictive. We extend the 
definition of distance-cumulative weighted graph properties X by assuming that X{G) = oo for 
every disconnected graph G. A weighted graph property X is said to be efficiently calculated if 
there exists a polynomial time algorithm that computes X{G) for every graph G. 

Consider some graph G G G. The diameter of G, denoted diam(G), is defined as diam(G) = 
max{distG(u, f ) \ u,v £ V{G)}, where distG{u,v) is the distance between u and v in G, namely, 
the length of a shortest (with respect to the edge weights) path from u to v. For a designated 
vertex u G V{G), the radius of G, denoted rad(G), is defined as rad(G) = iiiax{distG{u,v) \ 
V G V{G)}. The weight of a minimum spanning tree of G, denoted MST(G), is defined as 
MST(G) — mm{Y2e^T '^{^) I r is a spanning tree of G}. It is easy to verify that diam(G), rad(G), 
and MST(G) are distance-cumulative weighted graph properties. Efficient algorithms that com- 
pute them are described in most textbooks on graph algorithms (e.g., [12]). Another distance- 
cumulative weighted graph property is the diameter of a min- diameter spanning tree, defined as 
min{diam(T) | T is a spanning tree of G}. The min-diameter spanning tree is also known to be 
efficiently calculated [19] . 

Two more efficiently calculated weighted graph properties that fit into the definition of distance- 
cumulative properties by slightly modifying requirement (R5) are the all pairs average distance and 
the single source average distance with respect to some designated vertex. Our techniques can be 
adjusted to handle such a modification, however, this is beyond the scope of the current version of 

^ Graph G is said to be connected if it admits a path between u and v for every two vertices u,v £ V{G). 
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the paper. 



Our contribution. Let X be some distance-cumulative weighted graph property and assume 
that it is efficiently calculated. Given some connected RW graph G, X{G) is a random variable — 
denote it by X — and we are interested in approximating its fc*^ moment, i.e., E[X'^], for any fixed 
k. Specifically, we develop a fully polynomial time randomized approximation scheme (FPRAS) 
for the problem, namely, a randomized algorithm that runs in time poly(|G|, 1/e) for any choice 
of e > 0, where |G| stands for the number of bits required to encode G in a standard binary 
representation, and returns a (1 + e)-approximation of E[X'^] with probabilitj|§ at least 3/4. The 
following theorem is established in Sections [3l HI and [5j 

Theorem 1. The problem of computing the k^^ moment of X{G) on connected RW graphs G admits 
an FPRAS for every fixed k. 

In general, Theorem [1] is best possible. This is because exact solutions to the problems of 
computing E[diam(G)] and E[rad(G)] are ^^P-hard to obtain even when the input is restricted 
to identically distributed weighted graphs. Refer to Section [6] for a proof of this rather simple 
observation. 

To the best of our knowledge, our FPRAS yields the first provably polynomial time algorithm 
with guaranteed approximation ratio for any non-trivial statistical feature of a weighted graph 
property in randomly weighted graphs. Moreover, it seems that most related literature focuses 
on individual weighted graph properties and does not attempt to provide a framework for a more 
general theory of such properties; indeed, Snyder & Steele call for such a framework in their survey 
|32j . We hope that our technique which is suitable for all distance-cumulative weighted graph 
properties will be a significant step in that direction. 

Related work. The algorithmic aspects of randomly weighted graphs have been extensively 
studied since the early 60's (cf. Fulkerson [17J) mainly in the context of the shortest (s,t)-path, 
the longest (s,t)-path (a.k.a. the PERT problem), and maximum (s,t)-fiow. A comprehensive 
account of the various methods developed for the computation (and approximation) of the statistical 
features corresponding to these weighted graph properties is provided by Ball et al. [3] who also 
observe that an exact computation of the expected values of these weighted graph properties is 
#P-hard. Note that except for a few special cases (e.g., series-parallel networks with a specific type 
of edge distributions), none of the algorithms developed in that context is provably polynomial 
with guaranteed approximation ratio. 

Distance-cumulative weighted graph properties were also investigated under this setting. Hassin 
&: Zemel [20] prove that the diameter (and radius) of a complete n-vertex graph whose edge weights 

^ Using the median of means method, the success probabihty of an FPRAS can be increased to 1 — e for any choice 
of e > at the cost of increasing the run-time by an C'(log(l/e)) factor. 
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are uniformly and independently distributed in [0,1] is almost surely 0(log(n)/n). The hidden 
constants in this expression were resolved by Janson [23j who shows that the {s, i)-distance, radius 
with respect to s, and diameter of a complete n- vertex graph whose edge weights are uniformly 
and independently distributed in [0,1] converge in probability to Inn/n, 2 In n/n, and 3 In n/n, 
respectively. Frieze |14j shows that for every distribution function F with finite variance whose 
derivative at zero exists and satisfies F'{0) = D > 0, ii the edge weights in a complete n- vertex 
graph G are independently distributed according to F, then the weight of a minimum spanning 
tree of G converges in probability to (,{?>)/ D, where ("(3) = X^^^i 1/j'^- This is generalized by Steele 
|33j who shows that the assumption on the finite variance can be lifted. For the special case of F 
being the uniform distribution over [0, 1], Beveridge et al. [7] establish a variant of this bound for 
r-regular graphs. 

Asymptotic results for minimum spanning trees on n points uniformly and independently dis- 
tributed in the Euclidean unit ball are established by Bertsimas & van Rysin [6]. Kulkarni [28] and 
Alexopoulos & Jacobson [1] present algorithms that compute the distribution of MST(G) for graphs 
G whose edge weights obey exponential and discrete distributions, respectively. The run-times of 
their algorithms are not necessarily polynomial, though. A non-trivial upper bound on E[MST(G)] 
is established by Jain & Mamer |22j . 

A different, yet, related subject that admits a plethora of literature is average case analysis for 
graph algorithms (e.g. [151 [TTl There it is assumed that the edge weights in the input of 
some graph algorithm are drawn from a specified probability distribution and the goal is to analyze 
the expected run-time of the algorithm with respect to that distribution; refer to [16] for a survey. 
It is important to point out that upon invocation of the graph algorithm, the actual edge weights 
are determined, and in particular, known to the algorithm, in contrast to the RW graphs setting 
in which the challenge is to cope with the uncertainty in the edge weights. 

Techniques. We now provide an informal overview of the construction of an FPRAS for E[A'(G)], 
where X is some efficiently calculated distance-cumulative weighted graph property and G is an 
RW graph; approximating higher moments is very similar. Note that if the variance of X{G) is 
at most some polynomial times E[A'(G)]^ (that is, the critical ratio is polynomial), then E[Af(G)] 
can be approximated by means of sampling (a.k.a. Monte Carlo method). However, there exist 
some simple examples (cf. Appendix [X]) in which the variance is too large, and therefore a different 
approach must be sought. 

Our first step is to employ a simple reduction that allows us to focus on BCW graphs rather 
than arbitrary RW graphs. This reduction is presented in Section [SJ So, in what follows our goal 
is to approximate E[^Y(G)] to within a multiplicative error of 1 + 0(e), where G is a BCW graph. 

The desired approximation would have been straightforward to obtain if we could have efficiently 
approximated P(Af(G) > x) for arbitrary choices of real x > 0. It turns out that the case x = is 
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well known: approximating ¥{X{G) > 0) is equivalent to approximating the all-terminal network 
reliability problem (cf. |3j) in the graph obtained from G by removing all edges whose light 
phases are strictly positive. The FPRAS developed by Karger in [21] for the all-terminal network 
reliability problem is used to obtain this approximation. Karger's technique is based on identifying 
a collection C of polynomially many (2-way) cuts in G^ such that the probability that all edges of 
some cut not in C take on their heavy phases is small. (Of course, Karger uses the language of 
network reliability, rather than that of BCW graphs.) 

Unfortunately, approximating ¥(X{G) > x) to within a multiplicative error of 1 -1- 0(e) for an 
arbitrary choice of real x > seems to be a challenging task. In fact, when X = MST (namely, 
we are required to approximate the expected weight of a minimum spanning tree), an efficient 
implementation of this task would yield an FPRAS for the Tutte polynomial TG{x,y) of arbitrary 
graphs G for every x,y > 1 (refer to Bollobas [Sj for a comprehensive treatment of the Tutte 
polynomial and its many applications); whether or not such an FPRAS exists is an important open 
question piTS]. 

Instead, we use a careful "sliding window" argument to show that the desired approximation of 
E[Af(G)] can be efficiently obtained by repeatedly calling a procedure called Procedure Estimate. 
Given a positive real x, Procedure Estimate approximates ¥{X{G) > x) to within an additive error 
of 0(e) • ¥{X{G) > 0). In other words, we show that it suffices to produce a weaker approximation 
of ¥{X{G) > x); the quality of this weaker approximation is determined by ¥{X{G) > 0). The 
"sliding window" argument, presented in Section [3l is based on iteratively resetting all edge phases 
in G smaller than some threshold for a carefully chosen sequence of thresholds. 

Procedure Estimate itself is presented in Section HI The main trick there relies on formulating 
the real valued random variable Y that maps each instance I of the probability space defined by 
E{G) - E{G°) to F{X{G) > x \ I). Since E[Y] = F{X{G) > x), it is sufficient to approximate 
E[y]. We would have wanted to do so by sampling instances Z of the probability space defined by 
E{G) — E{G^) and then computing ¥{X{G) > x \ X). Sampling instances of the probability space 
defined by E{G) — E{G^) is a straightforward task. The problem is that given such an instance X, 
it is not clear how to efficiently compute P(A'(G) > x\X). 

To tackle this obstacle, we revisit the cut collection C and for each instance X, identify those 
cuts in C that conditioned on X, imply X{G) > x. For that to work, we must extend Karger's 
construction of C to r-way cuts for all r > 2. This extension builds upon the recent bound of 
Berend &: Tassa on the Bell number [5]. We then employ the method of Karp, Luby, and Madras 

[27] for probabilistic DNF satisfiability to approximate the probability that at least one of these 
cuts is induced by the edges in E{G^) that take on their heavy phase. 
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2 Preliminaries 



Randomly weighted graphs. Throughout we consider some distance-cumulative weighted 
graph property X and some n-vertex connected RW graph G. Let X denote the random vari- 
able that takes on X{G). By requirement (Rl), we may assume without loss of generality 
that the edge phases of G are scaled so that the smallest non-zero phase is exactly 1. Con- 
sequently requirement (R5) implies that X is either 0, or it is bounded from below by 1. 
On the other hand, requirements (R2) and (R4) guarantee that X is bounded from above by 
Xmax = ■ maxjl^g | e G E{G) and 1 < i < m(e)}. 

Suppose that G is a BCW graph. In what follows we assume that < pe < 1 for every edge 
e G E(G) (this assumption is clearly without loss of generality as the phases of an edge are not 
required to be disjoint). Let F C E(G) be some subset of the edges. Each edge e £ F takes on 
its heavy phase with probability pe and on its light phase with probability 1 — pe] this defines a 
probability space. It will be convenient to view an instance X of this probability space as a Boolean 
function I : F ^ {H , L} , where 



At the risk of abusing notation, we may sometime write F when we actually refer to the probability 
space it defines; our intentions will be clear from the context. 

Cuts and compact cuts. Consider some connected graph G. An r-way cut C of G is a partition 
of V{G) to r pairwise disjoint subsets, that is, G = {Ui, . . . ,Ur}, where Ui<j<r ^« ~ ^i^) and 
Ui n C/j = for every i ^ j- The subsets Ui, . . . ,Ur are referred to as the clusters of C. A cut 
refer^ to an r-way cut for any r > 2. 

Consider some r-way cut G = {Ui, . . . , Ur} of G. We say that an edge e G E{G) crosses G if 
e G UiX Uj for some i ^ j- The set of edges crossing G is denoted by E{G). The cardinality \E{G)\ 
is referred to as the size of C; if the edges of G are assigned with positive costs c : E{G) — )• M>o, 
then the sum "^^,^^(0) ''(^) referred to as the cost of G. The cut G is called compact if G{Ui) is 
connected for every 1 < i < r. Note that every r-way cut is a compact r'-way cut for some r' > r. 
A min cut (respectively, min cost cut) is a cut of minimum size (resp., cost). It is easy to verify 
that a min cut (resp., min cost cut) must be a compact 2- way cut. 

Consider some subset F C E{G) and a compact cut G of G. We say that F induces the cut G 
if the connected components of the graph obtained from G by removing the edges in F agree with 
the clusters of G. In particular, F must be a superset of E{G); it may contain additional edges as 
long as the removal of these edges does not disconnect any cluster of G. 

* In some literature, a cut refers to a 2-way cut, while an r-way cut for r > 2 is called a multiway cut. We do not 
make this distinction. 
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Analogy to the all-terminal network reliability problem. Consider some BCW graph G 
and some distance-cumulative weighted graph property X . Suppose that we wish to approximate 
the probability that X{G) > 0. Let = {e G E{G) \ = 0} and let G° be the restriction of G to 
the edges in . Requirement (R6) implies that the event X{G) > depends only on the probability 
space E^; specifically, X{G) > if and only if X{G^) > 0. Moreover, if G^ is disconnected, then 
X{G^) = oo with probability 1, thus we subsequently assume that G^ is connected. By employing 
requirement (R6) once more, we conclude that X[G^) > if and only if the edges that take on 
their heavy (positive) phases under the probability space E^ induce a cut on G^. This leads us 
to an analogy between the problem of approximating F{X{G) > 0) and the all-terminal network 
reliability (ATNR) problem. 

The input to ATNR is a connected undirected graph H in which each edge e fails (i.e., removed) 
with some specified probability. The goal is to compute the probability, referred to as the failure 
probability of H and denoted FAIL(//), that H becomes disconnected following such an edge failure 
experiment. ATNR is known to be jjP-complete [341 [3] and Karger develops an FPRAS for it |24] . 
Since H becomes disconnected if and only if the failing edges induce a cut on it, we conclude that 
W{X{G^) > 0) = FAIL(G''), where the latter is defined over an instance of ATNR in which each 
edge e G E{G^) fails with probability pe (the probability that e takes on its heavy phase in the 
BCW graph framework). Consequently, approximating V[X{G) > 0) can be performed by a direct 
application of Karger's algorithm. 

Monte Carlo method and approximators. Consider some probability space with sample 
space n and let A : (7 — t- M be a real valued random variable over this probability space. Suppose 
that the expectation of X is defined and denote it by fi. Let Ai, . . . , A„ be n independent samples 
of A and fix A = X^"=i Aj/n. Evaluating by A is referred to as the Monte Carlo method (cf. 
|26j). Let e and e be some positive reals. The following two theorems are direct consequences of 
Chernoff's inequality [10| (Theorem 12. ip and Hoeffding's inequality [21] (Theorem 12. 2p . 
Theorem 2.1. If X is an indicator random variable (namely, X G {0, 1} then taking n > 
41n(2/e)/u/e^ samples guarantees that P(|A — /i| > e) < e. 

Theorem 2.2. // A is almost surely in the interval [a,b], where b — a = p, then taking n > 
ln(2/e)p^/(2e^) samples guarantees that P(|A — ^| > e) < e. 

This leads us to notion of approximators. Consider some non-negative real value v that we 
would like to approximate with the real v' . Then v' is said to be an {e,e)- approximator of v if it 
satisfies the inequality |u — v'| < e with probability at least 1 — e, where the probability is taken 
over the randomness used to generate v'. Under this terminology. Theorems 12.11 and 12.21 provide 
sufficient conditions to guarantee that A is an (e, e)-approximator of /i. 

Proposition 2.3. If vi is an {ei,ei)- approximator of vq and V2 is an {€2, €^2) -approximator of vi, 
then V2 is an (ei + £2,61 + 6^2) -approximator of vq. 

Proposition 2.4. If v'- is an {e,i)- approximator of Vi for every 1 < i < n, then X^ILi'^^i/"' ^'^ ^''^ 
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(e, n • e)- approximator of "^^^i Vi/n. 

3 An FPRAS for BCW graphs 

In this section we consider some n-vertex BCW graph G and some smah performance parameter 
e > 0; our goal is to approximate E[X'^] to within a multipHcative error of l + 0(e). Here, we restrict 
our attention to the case k = 1, that is, we approximate E[X]. Extending our result to larger (yet 
fixed) values of k is mainly a matter of notation and we omit it from this version of the paper. The 
approximation presented in this section builds upon the more sophisticated Procedure Estimate 
which is presented in Sectional 

Theorem 3.1. There exists a randomized algorithm that with probability at least 3/4, approximates 
K[X] to within a multiplicative error of 1 + 0(e) in time poly(|G|, 1/e). 

Let be the smallest integer such that Xmax < (1 + e)^. (Note that A'' is proportional to 
log(2;inax)/e = poly ( | G| , 1/e)) . Clearly, we have < X < {1 + e)^. Fix vr^ = F{X > (1 + e)^) for 
every < i < N . Towards the approximation of IE[X], we first define 

N 

\ ^/'l I — 1 TO /'/'I I 



A = 2^(1 + e)*-i • P ((1 + e)*-i < X < (1 + e)*) 
1=1 
N 

= ^(1 + e)-i • [vr,_i - vr,] 
1=1 

N-l 

= vro + e • (1 + ef-' • vr, - (1 + ef~' ■ tt^ 

i=l 
N-1 

= vro + ^ e • (1 + ef-' ■ , (1) 

i=l 

where ([1]) is due to the fact that tt^v = 0. It is easy to verify that 

E[X]/(1 + e) < ^ < E[X] , (2) 

so our next goal is to approximate A. Note that ([1]) enables the computation of a (1 + e)- 
approximation of E[X] based on (1 + e)-approximations of P(X > x) for sufficiently many values of 
X. Unfortunately, we do not know how to obtain such an approximation directly and we are forced 
to apply some modifications to G. 



The shrunk graphs. Fix k = log^^^^ ^!Lil±£)^ 



For i = 0, 1, . . . , iV - 1, let Gi be the BCW 
graph obtained from G by setting Wt ^ for every edge e G E{G) and phase ip G {H, L} such 
that We < (1 + ey~'^. We refer to a phase that was set to in this process as a shrunk phase; 
the graphs Go, . . . , Gn~i are called the shrunk graphs. (Clearly, if the heavy phase of some edge 
is shrunk, then so is the light phase.) 
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Let Xi be the random variable that takes on X{Gi). Requirement (R2) guarantees that Xq > 
Xi > • • • > Xjv_i. The assumption that the minimum positive phase is scaled to 1 implies that 
G = Gq = Gi = ■ ■ ■ = Gk, hence X = Xq = Xi = • • • = X^. If i > k, then Xi may be smaller than 
X, however it is not much smaller as depicted in the following proposition. 
Proposition 3.2. //X > (1 + e)* for some <i < N - 1, then X/{1 + e) < Xi<X. 

Proof. Since phase We shrinks in Gi only if We < (1 + e)*~'^ < (1 + e)* • , requirements (R2) 

and (R4) guarantee that 

The assertion follows. □ 
Fix it'- = F{Xi > (1 + e)^) for every < i < iV - 1 and define 

N-l 



1=1 

Clearly, vr^ < tTj for every < z < — 1, thus A' < A. On the other hand, for every 1 < i < N — 1, 
we have 

ttU = > (1 + e)*"') > > (1 + ey-^) > F{X > (1 + e)*) = tt^ , 

where the last inequality is due to Proposition 13.21 Since, ttq = ttq, we get 

e) • TTj 



, N-l 

^ ' 1=1 

N-l 

i=2 
N-2 

= vr^+ J^e-a + er^-^^ < A', 



therefore 



A/{1 + e)<A'<A. (3) 
So, our next goal is to approximate A' . 
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Relying on local approximators. Consider some n-vertex BCW graph H and denote the 
random variable that talces on X{H) by Xh- Let 5, 5 > Q he some performance parameters. Recall 
that Karger's FPRAS for ATNR can be used to generate a real 8^^ that serves a.s a. {5 ■ ¥{Xh > 
0), (5)-approximator of P(X/f > 0) in time poly(n, log(l/5)) (see Section [2]). In Section [4] we 
present Procedure Estimate that given some x > 0, runs in time poly(n, 1/5, log(l/5)) and returns 
a real that serves as a (5 • f{XH > 0), (5)-approximator of F{X{j > x). 

Set the performance parameters 5 ^ f n/i+ \ ) and 5 ^ 1/{8N). For i = 0, 1, . . . , A'' — 1, we 



invoke Karger's FPRAS and our Procedure Estimate on H Gi with x {1 + e)' to produce the 
local approximators — )• and — )• £^^. We then define 



^^0^° if i = 0; ;v-i 

max{fp, f^>Oj if < i < A - and ^" = < + ^ e • (1 + e)^"^ • < . 

£r'' if A - K < i < A - 1 



How well does A" approximates A'? In attempt to answer this question, we first establish the 
following lemma. 

Lemma 3.3. With probability at least 3/4, we have 

(a) IvT' — tt'/\ < Stt'q for all < i < k; and 

(b) - <| < 5^[_^ for allK<i<N- 1. 

Proof. By the choice of 5 = 1/(8 A), we may use a union bound argument and conclude that the 
inequalities 

\F{Xh>0)-£>°\<5-¥{Xh>0) (4) 
\F{Xh>x)-£^''\<5-¥{Xh>0) (5) 

hold (simultaneously) for all A invocations of Karger's FPRAS and Procedure Estimate with 
probability at least 3/4; the remainder of the proof is conditioned on that event. By the definition 
of the shrunk graphs, we have 

F{Xi > 0) = P(Ao > 1) = TTo for every 0<i<K (6) 

and 

P(A, > 0) = F{Xi > (1 + e)*^*") < P(Ai_« > (1 + e)^"^) = tt-^^ for every k < i < A - 1 . (7) 

Combining ([6]) and (jl]) implies that Ivtq — ttqI < 5tt'q as required for i = 0. 
Next, ([5]) guarantees that 

|P(A, > (1 + e)') -£f''\<5- F{Xi > 0) (8) 
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for every < i < — 1, while (j3|) guarantees that 

£>,l <P(X,+« > 0) + 5 • P(X,+« > 0) 

<P(X, > (1 + e)^) + 5 ■ P(X,+« > 0) (9) 
<¥{X, > (1 + e)^) + 6 ■ ¥{X, > 0) (10) 

for every < i < N — k, where ([9]) is due to d?]) and (jlOp follows from the definition of the shrunk 
graphs. By combining ([8]) and (jlOp . it follows that for every < i < — k, we have 

> (1 + eY) - 6 ■ FiXi > 0) < max{fp,£:,>Oj < P(Xi > (1 + e)^) + 6 ■ P(Xi > 0) , 

or in other words, 

|P(X, > (1 + e)^) - max{fp,£:>Oj| < 6 ■ ¥{Xi > 0) . 

This yields the desired — vrf | < Sttq for every < i < k due to ([6]); and Ivr^ — 7r"| < 6ttI_^ for 
every k < i < N — k due to ([7]) . 

It remains to show that |vr^ — vr^'l < ^'^'i-K every N — k < i < N — 1. This is a direct 
consequence of ([8]) and ([7]). □ 



We are now ready to complete the analysis. Lemma 13.31 guarantees that 



\A' - A"\ 



N-l 



i=l 
N-l 



< K-7r'^\ + Y.e-{l + ey-'-\7r',-4\ 



i=l 



< Stt' 



l + 5^e-(l + eri 



i=l 



N-l 



1 + e 



(l + e)'^-l 



i=K+l 

N-l 



i=K+l 



<5- (1 + e)'^ 



N-K-l 



1=1 



with probability at least 3/4. By the choice of 6 = ^^^^j^j:^^ , we conclude that \A' — A"\ < eA'. 
Theorem 13.11 follows by combining the last inequality with ([2]) and ^ . 
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4 Procedure Estimate 



In this section we present and analyze Procedure Estimate. Let G be some n-vertex BCW graph 
and denote the random variable that takes on X{G) by X. Consider some positive real x and perfor- 
mance parameters e, e > 0. Given G and x, Procedure Estimate runs in time poly(n, 1/e, log(l/e)) 
and outputs an (e • V[X > 0), e)-approximator of ¥{X > x). 

Recall that requirement (R3) implies that edges of zero weight can be contracted without 
affecting the value of X. This means that we may contract every edge e G E{G) such that 
= W^ = without affecting F{X > 0) and F{X > x). Indeed, in what follows we assume that 
> for ah edges e e E{G). 

Fix E = E{G). It will be convenient to partition the edges in E according to their light phases 
to E" = {e £ E \ = 0} and to E - E" = {e e E \ > 0}. Let G° be the restriction of 
G to the edges in E^. Recall that the event X = X{G) > occurs if and only if X{G^) > (see 
Section [2]) ; denote the probability of this event by Pq . 

Fix c = (5 + vT7)/2 4.56. If Pq is sufficiently large, specifically, at least then the desired 
approximation can be obtained by a direct Monte Carlo method. Indeed, Theorem 12.11 guarantees 
that a Monte Carlo method with 0(log(l/e)n^/e^) = poly(n, 1/e, log(l/e)) trials suffices to generate 
an (e- Pq, e)-approximator of P(X > x). (The random variable for which we apply the Monte Carlo 
method is simply the indicator of the event X > x.) This applies in particular to the case where G^ 
is disconnected which means that Pq = 1- Therefore in what follows we may assume that Pq < iT' '^ 
and in particular, that G^ is connected. Note that if Pq is extremely small (e.g., exponentially 
small in n), then the above Monte Carlo method requires too many samples in order to obtain an 
{e ■ Pq, e)-approximator of P(X > x). 

Dealing with small Pq. How do we efficiently generate an {e ■ Pq, e)-approximator of ¥{X > x) 
when Pq is small? For that purpose we introduce the real valued random variable Y which is 
defined over the probability space E — E^ by mapping each instance I : E — E^ — t- {H,L} to 
P(X > X \ I), namely, I is mapped to the probability that X{G) is at least x conditioned on 
I. This can be viewed as decomposing the probability space E into the Cartesian product of the 
probability spaces E — E^, from which I is chosen, and E^, over which F{X > a: | X) is defined. A 
crucial observation here is that 

E[Y] = F{I) ■ F{X >x\I) 

I:E-E0-^{H,L} 

F{X >xAl) 

I:E-EO^{H,L} 
= F{X > x) , 
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hence our goal is to provide a good approximation for E[y]. Another important observation is that 

¥{X >x\I) < F{X > I X) = F{X >0) = Pq 

for every instance I : E — — )• {H,L} (recall that the event X > d does not depend on the 
probability space E — E^), therefore Y £ [0, Pq] with probability 1. 

Fix k = 21n(4/e)/e^ and repeat the following process for j = 1, . . . ,k. Choose some instance 
Ij : E - E° ^ {H, L} with probability P(X,) (th is can be easily generated by randomly choosing 
the phase of each edge in E — E^ independently of all other edges) and let Yj = ¥{X > x \ Ij); 
in other words, Yj is a random sample of Y. Unfortunately, we do not know how to efficiently 
compute the exact value of Yj for a given instance Ij. Instead, we will generate an approximate 
sample Yj which is an (e • i-b/2, e/(2A;))-approximator of Yj. 

We will soon explain how the YjS are generated, but first let us explain how they are employed 
to obtain the desired (e • e)-a'Pproximator of K[Y]. Let Y = X])=i^/^ ^.nd Y' = Yl^j^iY- /k. 
Theorem 12.21 guarantees that Y is an (e • Po/2, e/2)-approximator of E[y]. By Proposition 12.41 we 
conclude that Y' is an (e • Po/2, e/2)-approximator of Y . Therefore Proposition 12.31 implies that Y' 
is an (e • Pq) e)-approximator of E[y] = P(X > x) as required. 

Generating the approximate samples. It remains to present the process through which the 
approximate samples Y- are generated (recall that each approximate sample should be an (e • 
Po/2, e/(2fc))-approximator of ¥{X > x \ T) for some given instance I : E — E^ ^ {H,L}). The 
technique we use for this process is an extension of Karger's technique [Ml- In order to simplify 
the description of this process, we first assume that there exists some real p such that pe = p for all 
edges e £ E^ (all these edges have zero light phase, however their heavy phases may vary). This 
assumption is removed later on. 

Given some compact cut C of G*^, let 'H(C) denote the event that all edges in E{C) C E^ take 
on their heavy (positive) phases. Let C be some min cut of and let x = l-^(C)l be its size. Since 
^("^((7)) = p-^, the assumption that Pq < n~'^ implies that p^ < n"'^. The following two theorems 
are established in [23] for the case of 2-way cuts. Building upon the recent bound of Berend & 
Tassa on the Bell number [5j, we extend them to (compact) r-way cuts for all r > 2 simultaneously. 
Theorem 4.1. For every real a > 1, there are less than 13n^° compact cuts of size at most ax in 
. Moreover, these cuts can he enumerated in expected time 0(n^"). 

Proof. The theorem is established by presenting a random process that generates each compact cut 
of size at most ax in G^ with probability greater than Observe first that if an r-way cut 

C satisfies |£'(C)| < ax, then r must be at most 2a as otherwise there exists some cluster U of 
G with less than x edges crossing between U and V{G^) — U, in contradiction to the assumption 
that X is the size of a min cut of G^. 
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Fix k = [2a] . Our random process first performs random edge contractions in until k 
vertices vi,...,Vk remain in tlie grapli (cf. Section 2.2.1 in |24j): eacli vertex Vi corresponds to 
some subset Vi C V{G^) so that {Vi, . . . , Vk} is a partition of V(G^) (the subgraph induced on G^ 
by Vi is connected). We then take P = {Wi, . . . , Wi}, 1 < £ < A;, to be a partition of {vi, . . . , Vk} 
chosen uniformly at random out of the Bj^ possible partitions of {vi, . . . ,Vk}, where Bf^ is the k^^ 
Bell number. The cut C = {Ui, . . . ,Uf} generated by our random process is defined by setting 
Uj = [j{Vi I Vi G Wj} for j = 1, . . . , i. It is important to note that C is not necessarily a compact 
cut, however, if C is any compact cut of size at most ax in G^, then it can be generated by our 
random process. Our goal in the remainder of this proof is to show that C is indeed generated with 
probability greater than 

Karger [24J shows that the probability that none of the edges crossing C is contracted during 
the random edge contractions is at least 



\2a) 

where generalized binomial coefficient^fl are used when 2a is not an integer. It remains to prove 
that i^B-^ > Indeed, 



(2J 0-1 ^ ( k\ 2a\ „ 1 



Q"J^^ " \2a) \en) 







> 


(- 




\2a 






> 


(en. 






> 





2a 



0.792fc 

k 



1 /ln(A; + l 



k V 0.792e 

where inequality [11] is due to Berend & Tassa [5]. The assertion follows as ^ (^^^^^^^ > 1^ when 
k>2. □ 

Theorem 4.2. For every real a > 1, the probability that there exists some compact cut C of size 
at least ax in G*^ such that TiiC) occurs is 0{n~^'^), where rj is defined by fixing = n~^'^~^^\ 

Proof. Let Ci, . . . , Cj be the compact cuts of size at least ax and for each 1 < i < t, let Xi = \l^{Ci)\. 
We assume without loss of generality that ax < Xi ^ " " " ^ Xt- Denote pi = p^^ = P('H(Cj)) and 
consider some real /3 > 1. By Theorem l4.lt there are less than ISn^'^ compact cuts of size at most 
j3x- It follows that Xisn^/s niust be greater than f3x- 



^ Generalized binomial coefficients are a generalization of the standard binomial coefficients Q to non-integral x 
and y. This generalization is based on replacing the factorial in the standard definition with the Gamma function. 
Many of the identities and bounds that hold for the standard binomial coefficients also hold in the generalized case, 
including the bounds (^|)^ < < (^^y 
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We shall bound the probability that T-L{Ci) occurs for some (at least one) 1 < i < t by bounding 
the sum 

Y!i=iPi- The first t' = i3n2"+i/inW 

terms are bounded simply by observing that 

t' 

^Pi< 13n2°+Vln(n) . pax ^ . ^2a . ^-a(2+r,) ^ -^3^ . ^-arj _ 

i=l 

Thus it remains to bound the remaining t — t' terms. 

Given some /3 > a, we write s = I3n'^l^+^/ ^^("'^ and conclude that Xs > (^13+ 2 in{n) ) ^ ~ 

ln(s)— ln(13) rpi r 

2 Inn — ■ /t- J-heretore 



Summing over all i > t' , we get 



< 13l+V2.^^-{l+V2) 
i>t' s>t' 

POO 



V 13n2a + l/ln(n) 

= 13l+''/2 . (^/2) • 13-(''/2) . ^-av . g-r,/2 

< 13 • n-"'' . 

The assertion follows. □ 



Notice that the compact cuts addressed in Theorems 14. II and 14. 21 may have an arbitrary number 
of clusters, but their size is compared to ax, where x is the size of the smallest 2-way cut in G^. 
This point is crucial for the validity of the arguments. 

Write = n~(2+''). We must have rj>c — 2asp^< n^'^. Fix a = !"("■) g^j^^j jg|- 

C be the collection of all compact cuts of size at most ax in . By Theorem 14.11 C consists of 
0(71.2°) _ (^(T^c-i^g^ compact cuts that can be enumerated in expected time 0{n^~^ / e). 

Given some sub-collection i3 C C, let 

^{B) = ^[\l n{c)\ 

\ceB ' 

be the probability that all crossing edges of some (at least one) cut in B take on their heavy phases. 
Theorem 14.21 guarantees that < Pq — V'(C) < 771""'' for some universal constant 7. The choice 
of c = (5 + ^/Y1)/2 and of a = ^ ^+^°(y^)/ ^"(") ^^^^ assumption that 77 > c — 2 ensure that 
r^ji-otr] ^ en~(2+'?)y'4 = as long as (47)^ < (l/e)^~^, which yields the following corollary. 

Corollary 4.3. The probability that there exists some compact cut C ^ C such that 'H{C) occurs 
is at most ep^/A < e ■ Po/4. 
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Consider some instance X : E — — t- {H,L}. Our goal is to efficiently generate an (e • 
Po/2, e/(2A;))-approximator of P(X > x | X). For a given compact r-way cut C = {Ui, . . . , Ur} of 
G^, we construct the graph Gc,i as follows. The vertex set of Gc,x is V{Gc,i) = {ui, . . . , Uj-}. For 
every edge e ^ E^ with one endpoint in the cluster Ui and the other in the cluster Uj, i 7^ j, we 
add an edge (u^, Uj) to E{Gc,x) whose weight is . In addition, for every edge e ^ E — E^ with 
one endpoint in the cluster Ui and the other in the cluster Uj, i ^ j, we add an edge {ui,Uj) to 
E{Gc,x) whose weight is We^^\ The following observation is due to requirement (R3). 
Observation 4.4. Conditioned on the instance I : E — E^ — t- {H,L}, and on the event that the 
set of edges in Eq that take on their heavy phases induces the compact r-way cut C on G^ , we have 
X{Gc,x) = X. 

Let Bx be the collection of all compact cuts G of G^ such that X{Gc,x) ^ x. Observation 14.41 
implies that ¥{X > x \ X) = iIj{Bx)- By Corollary 14.31 we know that iIj{Bx — C) < e-Po/4, and hence 
V'(%) - e • Po/4 < tl^iBx n C) < ipiBx). Consequently, it suffices to generate an (e • Po/4, e/{2k))- 
approximator of iIj{Bx n C). 

Probabilistic DNF satisfiability. The approximation of 'il){Bx^C) is performed by the method 
of Karp, Luby, and Madras |26|, I27j for approximating the probability that a formula in disjunctive 
normal form (DNF) is satisfied. Given some DNF formula (j), and given the probability qi that 
Xi is assigned to true for each variable Xi (independently of all other variables), the method of 
Karp et al. generates a {5 ■ g(0), (^)-approximator of the probability q{(j)) that is satisfied in time 
0(101 log(l/(5)/(5^), where |(/)| stands for the size of the formula (number of literals). 

Cast into that framework, the event VceBinc ^(^) encoded as a DNF formula whose vari- 
ables correspond to whether or not the edges in E^ take on their heavy phases and whose clauses 
correspond to the cuts in Bxf^C. Such a DNF formula has |;SinC| < \C\ = 0(n'^~^/e) clauses, each 
with at most n literals. Therefore an (e-'i/'(;SxnC)/4, e/(2A;))-approximator tl^{Bxr]C), which also 
serves as an (e • Po/4, e/(2A;))-approximator of ip{Bx n C) since tp{Bx n C) < Pq, can be generated in 
time 0{log{k/e)n'^/e^) = poly(n, 1/e, log(l/e)). 

Varying heavy phase probabilities. Recall that in attempt to simplify the description of the 
process that generates approximate samples of the random variable Y, we assumed that Pe = P for 
all edges e G E^. We now turn to lift this assumption. The technique we use here is essentially 
identical to that used by Karger [25j for a similar purpose; we describe it for completeness. 

The BCW graph G^ with varying heavy phase probabilities < pe < 1 is transformed into 
a BCW graph H, V{H) = V{G^), with uniform heavy phase probabilities p = 1 — 6 for some 
sufficiently small 9 > 0. For each edge e = {u,v) £ E^ with heavy phase probability < < 1, 
we introduce a bundle of k^ = \n{l/pe)/0 parallel {u,v) edges in H with the same heavy and light 
phases as e (recall that the light phase of e is zero) . The probability that all ke edges in this bundle 
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take on their heavy phase is (1 — which converges to pe as 6 ^ 0. By requirement (R4), 

it is sufficient to generate approximate samples for the random variable Y with respect to the graph 
H in the limit as — )■ 0. The technique we introduced earlier in this section is suitable for that as 
H has uniform heavy phase probabilities. In particular it is sufficient to enumerate all the small 
compact cuts C of H, identify those inducing X{Gc,x) > x for a given instance X : E-E^ {H , L}, 
and then approximate the probability that all crossing edges of at least one of them take on their 
heavy phases. 

Note that changing the parameter 9 scales the size of cuts in H without changing their relative 
sizes. We construct a graph H' , V{H') = V{G^), with positive costs on the edges by assigning cost 
ln(l/pe) to each edge e € E^ (the phases of the edges in E^ are ignored in the context of H'). The 
small cost cuts in H' correspond to the small sized cuts in H; they can be enumerated by known 
techniques. 

Given the small cuts in H that induce X(Gc,i) > x, we have to approximate the probability, as 
— )• 0, that all crossing edges of at least one of them take on their heavy phases. We already argued 
that this is exactly the probability that all crossing edges of at least one of the corresponding cuts 
in take on their heavy phases. Approximating this probability is done as before by constructing 
the appropriate DNF formula and employing the method of Karp et al. \26\ 127] . 



5 Transforming RW graphs into BCW graphs 

So far, we have developed an approximation for E[A:'(G)'^], where G is a BCW graph. In this section 
we are interested in extending our algorithm to the (general) case of RW graphs. This extension 
relies on an efficient transformation of any RW graph G into a BCW graph G' such that the random 
variable X(G') is stochastically equivalent to the random variable X(G) and \G'\ = 0(|G|). This 
is similar to a method presented by Mirchandani [30] (see also |3]). 

Let G be an arbitrary RW graph and consider some edge e G E{G). Recall that there exists some 
positive integer m(e) and some non-negative phases W}, . . . , W^"^^^ and probabilities p\, ■ ■ ■ ,pT'^^\ 
where Y^^l^ p\ — such that w{e) = Wl with probability p\ (independently of all other edges). 
We assume without loss of generality that if m{e) > 2, then the phases of e are distinct (identical 
phases can be merged into one) and ordered so that < • • • < W^"^^^ ■ 

The BCW graph G' is obtained by replacing every edge e = {u,v) G E{G) such that m{e) > 2 
with m(e) — 1 parallel edges ei, . . . , e^(g)„i = {u, v) £ E{G'), each having exactly two phases. The 
heavy phases of all new edges are set to be W^^'^\ The light phase of is set to be W^. = Wl for 
every 1 < i < m(e) — 1. The probabilities p{ei), . . . ,p(em(e)_i) are designed to guarantee that the 
random variable Me = min{?i;(ei), . . . , ^(6^(6)-!)} ™ G' is stochastically equivalent to the random 
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variable w{e) in G. This is achieved by setting 

pi 



P{ei) = 1 



Indeed, for every 1 < i < m(e), we have 



i-1 



p(Me = w^) = Hpie,) . (1 - p(e.)) = n ]_^r-/: ■ ti$^ = ' 

where the last equation holds by telescoping. Requirement (R4) implies that X{G') is stochastically 
equivalent to X{G). Theorem [T] now follows (at least when A; = 1) from Theorem 13.11 



6 Hardness 

In this section we prove that the problem of computing the expected diameter of an RW graph is 
#P-hard. The problem remains ^^P-hard even when restricted to identically distributed weighted 
graphs. Our line of arguments immediately implies that computing the radius of an identically 
distributed weighted graph with respect to a designated vertex is also #P-hard. 

Hardness is established by reduction from the most basic variant of the two terminal network 
reliability (TTNR) problem defined as follows. On input connected graph G and two vertices 
s, f e y{G)-, the goal is to compute the probability P that s and t become disconnected when each 
edge in E{G) is removed with probability 1/2 independently of all other edges. The #P-hardness 
of TTNR is established by Valiant [34]. Since the support of P consists of integer multiples of 2"™", 
where m = \E{G)\, we conclude that it is ^^P-hard to approximate P to within a one-sided additive 
error of e for any e < 2"*^. 

Given a graph G with two vertices s,t S ^(G) as input of TTNR, we construct an identically 
distributed weighted graph G' with parameters p = 1/2, = 1, and = 0. G' is obtained 
from G by adding a new edge e = (s, t) and augmenting the resulting graph with two simple paths, 
one connecting s to the new vertex s' and the other connecting t to the new vertex t' . Each new 
simple path consists oi k = Q{m) new vertices. The reduction is cast in the following lemma. 
Lemma 6.1. Let D he the random variable that takes on diam(G'). Then P < 2(E[D] — k) < 
P + 2-™. 

Proof. Let D' be the random variable that takes on distG'/(s', t'). By definition, we know that 
D' < D with probability 1. We shall take k to be sufficiently large so that Chernoff's inequality 
implies that P(distG/(s', s) < n V distcK*, < n) < 2'^'^+^'^ /n, where n = \ V{G)\. By the 
construction of G' , it follows that F{D > D') < 2-(™+iV'^- Since D — D' < n, we conclude that 
< E[D] -E[D'] < 2-(™+i). 
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By the linearity of expectation, we have 

E[D'] = E[distG'(s',s)] + E[dist(s, t)] + E[distG'(t, i')] = A: + E[dist(s, t)] , 
where the last term can be rewritten as 

E[dist(s, t)] = E[dist(s, t) \ w{e) = 1] ■ ¥{w{e) = 1) + E[dist(s, t) \w{e)=0]- F{w{e) = 0) 
= E[dist(s,t) I w{e) = l]/2 . 

The assertion is established by arguing that E[dist(s,t) | w{e) = 1] = P. Indeed, when w{e) = 1, 
then distG'{s,t) £ {0, 1}. The argument holds since the distG''(s,t) = instances (respectively, the 
distG'(s,t) = 1 instances) of the probability space E{G') correspond to the instances of TTNR in 
which s and t remain connected (resp., become disconnected). □ 

7 Conclusions 

We study the setting of graphs whose edge weights are independent random variables and show that 
for the wide family of distance-cumulative weighted graph properties, the problem of computing 
the k^^ moment admits an FPRAS. Computing the expectation (i.e., the first moment) for example 
is difficult when the variance is large, and hence too many samples are required for the Monte 
Carlo method in order to take into account low probability eventqj that may drastically affect the 
expectation. Our technique does not guarantee a (multiplicative) approximation for the k''^ central 
moment (and in particular, the variance) when this is close to zero, however, it does provide us 
with the ability to decide if the k^^ central moment is close to zero. 

There are still some fundamental weighted graph properties which are not distance-cumulative, 
and hence cannot be dealt with via our technique, such as the shortest (s, t)-path and the weight 
of a maximum matching. It is also natural to consider the directed analogue of randomly weighted 
graphs and in particular, various network flow problems. Another aspect that calls for further 
research, once encoding issues are resolved, is that of continuous distributions for the edge weights. 
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APPENDIX 



A High critical ratio 



Consider the identically distributed weighted graph G consisting of 2 vertices and m parallel edges 
connecting them, where each edge is of weight 1 with probability 1/2; and of weight 2^"* otherwise. 
Let X denote the random variable that takes on the diameter of G. It is easy to verify that 
E[X] 2"^, while Var[X] ^ 2^"*, so the critical ratio here is roughly 2"*. Indeed, a Monte Carlo 
method with significantly less than 2"* samples would probably estimate the expected diameter of 
G to be 1 which is an awful approximation. 
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